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' Abstract 

oo : 

We determine the path of the hght around a dielectric vortex described by the 
relativistic vortex flow of a perfect fluid. 

^ ■ 

o 

5 '. 1 Introduction 
o 

^ ■ Recently, Leonhardt and Piwnicki have studied the light propagation in a vortex flow of 
bJQ. non-dispersive dielectric medium in order to obtain a so-called optical black hole. However, 
Visser has shown, in a comment of another paper of Leonhardt and Piwnicki , that this 
vortex geometry has not the sense of black hole such that it is defined in general relativity. 



;h ' Moreover, we emphasize that their vortex fiow is simply an ad hoc generalization of the 
^ Newtonian case. Despite these points, the analysis of the light propagation in a dielectric 

vortex is interesting in principle. 

In this work, we take up again the question of the light propagation around a dielectric 
vortex in the case where the relativistic vortex flow is really determined within the frame- 
work of relativistic hydrodynamics. It is different from the vortex flow considered by the 
authors. Then, we establish the qualitative features of the trajectories of the light rays 
dependent on the values of the angular momentum of the light rays. 

The plan of the work is as follows. In section 2, we derive the relativistic vortex flow. 
For this case, we calculate in section 3 the effective metric for which the light rays deflne 
a null congruence of geodesies. In section 4, we obtain the equations determining the path 
of the light. Then, we discuss the qualitative features of the trajectories of the light rays 
in section 5. We add some concluding remarks in section 6. 



"E-mail: linet@cclfi.phys.univ- tours. fr 



2 Relativist ic vortex flow 



We are concerned with a perfect fluid having the energy- momentum tensor 

T"'' = {pc^ + p)u^'u'' + pri^"" (1) 

where p is the mass density, p the pressure, the four- velocity of the fluid, and rj^^ is the 
Minkowskian metric with signature ( — h ++). It is conserved d^jT^^ = 0. We assume the 
existence of a conserved mass density r, i.e. df^{ru^) = 0. We also define the relativistic 
specific enthalpy h = {pc^ +p)/r. 

We suppose that the motion of the fluid is isentropic. So, it is not necessary to consider 
the specific entropy and we postulate merely an equation of state h{p). From the equations 
of the fluid motion, it is easy to prove that dh/dp = 1/r. Moreover, by setting 

p^ = hu\ (2) 

we obtain the following equation 

u^drPa = —do-h. (3) 

We describe the Minkowskian spacetime in cylindrical coordinates {x^, p, (p, z) in which 
the Minkowskian metric has the form 

ds^ = -{dxy + dp" + p2rf<^2 ^ ^^2_ (4) 

The symmetries of the vortex flow are given by the three Killing vectors : do, d^, dz, denoted 
^(a) — 152,3), since we assume that the solutions to the equations of the fluid motion 
are stationary and cylindrically symmetric. In particular, we have ^(a)C?r^ = 0. Taking 
into account the fact that d^^(^a)u + dyi{a)^i = 0, we obtain from (H) that d„{^J^^^Pr) = 0. 
Thus, Cla)Pr are constants of the fluid motion. We limit ourselves to the case p^ = and 
we put 

Po = —E , /ip = , p^ = M. and Pz = k (5) 

where S {S > 0), Ai and k are arbitrary constants of the fluid motion. By using a Lorentz 
transformation along z, we may take k = 0. 

We now solve the equations of the fluid motion. From (^, we get h'^ = —prp'^- Hence 
by assumptions (|^), we find 

/ 

Kp) = f'--r (6) 

which is well defined for p > p^ where pc =\ M./E |. We notice that h{pc) = 0. By 
combining and (|^), we get the components of the four- velocity of the fluid 

M° = , ^ , = , m'^ = ± , and u' = for p> (7) 

- Pl/p' p'^^-pVp' 

in which we choose the sign + henceforth. Solution (|^) to the equations of the fluid motion 
gives a model of vortex flow which is only considered for p > p^, i.e. outside the core region 
of the vortex flow. This vortex geometry is usually used in relativistic hydrodynamics 
and in the relativistic theory of superfluid 
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3 Light propagation in a non-dispersive medium 

The electromagnetic properties of an isotrope non-dispersive dielectric medium are charac- 
terized by the permittivity e and the permeability /i. A relativistic model of such a medium 
has been found a long time ago P, |^ . The Maxwell equations are written with the aid of 
a tensor g^j^^, defined by 

9al3=Va(S- ( -1 + — ) UaUi3 (8) 



,2 _ ^,,„2 



where n is the refractive index, n = efic , and the four-velocity of the medium. At the 
approximation of geometrical optics in this medium, the vector of propagation of light is 
kfi = d^S where S is the phase and it satisfies the eikonal equation [^] 

r^Kkp = (9) 

where 

gag 

the inverse matrix of Qf^^, having from the expression 

g'"^ = 7]"^^ - [-1 + n^) u'^u^ . (10) 

Equation (^) for the propagation of light in this dielectric medium is equivalent to the 
propagation of light in the effective metric g^^. We therefore define the components k'^ of 
the vector of propagation by 

p = r^k, (11) 

and we take fc^ = k^. Since k^ = (9^5, equation (^ implies that the vector of propagation 
k^ defines a null congruence of geodesies for the effective metric g^y. We have thereby 

FV/,r = (12) 

where V is the metric connection associated to gfj,^. 

We suppose that the dielectric medium consists of a perfect fluid. In the case of the 
vortex flow given by (|^, we can calculate the components of the effective metric 

~ 1 - n'pl/p' . Pc (1 - l/n') . (1 - pg/nV) 

cjpp = 1 and g^z = 1 
and these of the inverse metric 



-00 _ ^ Hc/H -Oy _ Hc\ -^-r J ~w _ " f^c/ r /..x 

^ 1 - p2/p2 ' ^ p2 (1 _ ^2/^2) ' ^ ^2 (1 _ ^2/^2) ' ^ ) 

gPP = 1 and g'' = 1. 

Expressions ( |13D and (0) are well defined for p > pc, i.e. outside the core region of the 
vortex flow. 

The three vectors do, and dz define again three Killing vectors of the effective 
metric (jf^u, i-e. Va^{a)i3 + ^i3^{a)a = 0. Taking into account the geodesic equations ([T^), 



'-Pc/P' ^o<p Pc(-l+n2) l-n2p^/p2 



we find that ^{^^^^ku) = and so i^^^^k^ are constants of tlie liglit motion and tlius we 
put 

fco = —E k^ = L and k^ = C (15) 

wfiere E [E > 0, L and C are arbitrary constants, E being tlie energy and L the angular 
momentum of the hght rays. 



4 Equations of the path of the hght 

We parametrize the path of the light a;'^(A) with the affine parameter A so that we have 

obeying the geodesic equations (|T^ for p > pc. We limit ourselves to consider trajectories 
of the light rays in the plane orthogonal to the vortex flow. We can put ^(A) = having 
C = 0. We do not write down the geodesic equations ([T2| ) but we know that these are 
differential equations of second order. To determine uniquely the path of light, we must 
know at A = Ao the quantities 



However, there exists the constraint (P) on these initial data and consequently there is five 
independent initial data. 

We consider the characterization of the solutions to the geodesic equations (0) by the 
constants ([T5|). According to (0), we have the differential equations 



^ = -~9''E+~g'^L and ^=r^L-~g'-E. (17) 
aA cLA 

The unknown fonction p(A) is determined from equation (P|) which can be written 



dp 
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-g^^E^ + 2g°^EL - g^^L\ (18) 



In this case, the path of light is characterized by the five quantities 

Mo ' (<^)o ' (p)o ' E andL. 

Now we point out that the affine parameter A is only defined up to a constant factor, 
therefore we can put E = 1 hj a choice of the affine parameter. 

We now turn to analyse the differential equation (|TB]) for p > p^. By setting 

X = — , r = — and 7 = — , (19) 
Pc Pc L 
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equation ( ]T8| ) becomes 



dr 



X 



[x^ 



1)7^ 



n 



x^ 



7 



2(-l + n^ 



x^ 



-7 



n 



X^ / X^ 



(20) 



valid for x > 1 and 7 7^ 0. We can rewrite (^) in the form 

[(7-7i(x)) (7-72(3;))] 



dr 



(X2 - 1)72 



where 



n — X 



and 72(3;) 



n + z 



x{nx — 1) x(nx + 1) 

By virtue of ([iTl), the polar angle v^(t) is then determined by 

,2\ 1 



d(p 
dr 



1 - ^ I - - 1 + 

7 



f211 



(22) 



(23) 



when x(r) is known. 

In the case 7 = 1, we can easily find the general solution to equations (pT]) and (p3|) 



a;(r) = y^l + ?T,2(ro — r)^ and V9(r) = arctann(ro — r) + y^o- 



(24) 



In the case 7 7^ 1, expressions (^Tj) and (p3D have a divergence at x = 1 and therefore 
we have 



lim — 

\ dr 



00 and lim 

2:^1 \ dr 



00. 



Since the effective metric (|T3D is singular at x = 1, we cannot discuss the path of the light 
in terms of r at x = 1. In any way, the affine parameter A has no physical meaning since 
is an auxiliary metric. Also, it is convenient to express x in function of for x > 1 
and we obtain thereby 

( dx\ 
lim — = 0. 

\d({) ) 

This means that the trajectories of the light rays arriving at x = 1 are tangent at the 
circle x = 1. Then, we can match them with the similar trajectories outgoing at the same 
position. 

In consequence, we would directly determine x{ip) by the following differential equation 



' dx\ (x^ — l)(n^x^ — 1) 

d^J (1-^2)2 [7-70(0;)]^ 



[(7-7i(x))(7-72(x))] 



where 



x2 — 11? 

x^Cl — r? 



(25) 



(26) 



deduced from (PTI) and 



5 Qualitative features of the trajectories 



We are now in a position to discuss the qualitative features of the trajectories of the hght 
rays. Already for 7 = 1, we have the particular solution ( P^ which can be rewritten 

x(^) = ^l+(tan^^)^. (27) 

We see that the path of light described by trajectory (|27|) arrives at x = 1 from infinity 
then goes back at the infinity. 

A solution to equation (^) exists when its right member is positive. This fact depends 
on the value of 7 with respect to the roots 71(0;) and 72 (a;)- The turning points occur when 
dx/dif = 0. The derivative dx/dip changes of sign there except if the derivative of the 
right member with respect to 99 vanishes also at this point. In this latter case, we have a 
limiting circle. Moreover, as proved in the previous section, we see directly from equation 
(^) that the trajectories of the light rays passing at x = 1 are tangent at the circle a; = 1. 

In order to discuss graphically the form of the trajectories dependent on the values of 7, 
we must draw the curves 71 and 72. We note that 7i(l) = 1 and 72(1) = 1 and 7i(c>o) = 
and 72(c)o) = 0. The derivatives of 71 and 71 with respect to x are 

d'^i _ n{l - 2n + x^) 6/72 _ _n(l + 2n + x^) 

dx x^(nx — 1)^ dx x'^{nx + lY 

The derivative of the function 72 is always negative for x > 1 and so 72 is monotonically 
decresing from 1 to 0. On the contrary, the form of the curve of 71 depends on whether 
n < 1 or n > 1. 



5.1 Case n < 1 

The function 71 blows up for x = The derivative of the function 71 is always positive. 
For X > 1/n, we have 72(x) > 7i(x) and for 1 < x < 1/n, we have 72(x) < 7i(x). We 
remark that the component ^00 of the effective metric (0) vanishes for p = Pc/n and 
changes of sign. 

The solution to equation (|25| ) exists when 7 > 71 (x) or 7 < 71 (x) for x > 1/n and 
when 72 (x) < 7 < 71 (x) for 1 < x < 1/n since the factor in front of 7^ in (EBI) is presently 



negative. We see graphically on Fig. 1 that for all 7, 7 7^ 1, there is a turning point. We 
have quasi-hyperbolic trajectories. The case 7 = 1 is the particular solution (|27D already 
studied. 



5.2 Case n > 1 

The derivative of the function 71 vanishes at x = x/ with xi = n + ^n^ — \. We denote 
7i = which has the expression 
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7 



7 




Figure 2: Curves 71 (x) and 72(2;) for n = 1.2 



satisfying —1 < 71 (x;) < 0. We have 72(2;) > 71 (x). 

The solution to equation ( P^D exists when 7 > 72(0;) or 7 < 71 (x). We see graphically 
on Fig. 2 that when 7 > 1 or 7 < 7^, the light reaches tangentially to the circle x = 1 and 
goes back at the infinity. We point out the existence of a limiting circle at x = xi and thus 
for "J = '~fi the corresponding light ray spirals around this limiting circle. For 7; < 7 < 1, 
the form of the curves 71 and 72 shows that there are two turning points. We have either 
quasi-hyperbolic trajectories or quasi-elliptic trajectories, oscillating between x = 1 and x^ 
defined by 71 (x^) = 7. 



6 Conclusion 

In the present work, we have considered that the dielectric vortex is the relativistic vortex 
flow of a perfect fluid. The core region defined hy p < pc is not generally considered. 
Outside the core region, it is possible to calculate an effective metric for which the light rays 
define a null congruence of geodesies. By discussing the general features of the trajectories 
of the light rays, we have proved that the light rays remain external to the core region. 
The trajectories arriving at p = pc are tangent at the circle p = pc- However, we have 
merely an idealized model of dielectric vortex and, since the effective metric is singular at 
p = Pc we cannot determine the behaviour of the light within geometrical optics in the 
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neighbourhood of p = pc- 
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